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Abstract 

The effect of the cosmological constant on the curvature of light due to an isolated spherical 
mass is recalculated without using the lens equation and compared to a lensing cluster. 
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1 Introduction 



Until September 2007, the official opinion had it that the deflection of light by an iso- 
lated spherical mass is independent of the cosmological constant. This was corrected in 
the beautiful work by Rindler & Ishak [1]. One month later, Ishak, Rindler et al. [2] 
showed that this dependence is not academic by exhibiting several lensing clusters, that 
do constrain the cosmological constant within today's bounds from the standard model 
of cosmology 

Two consecutive papers by Lake [3] and Sereno jl] confirm the dependence. Sereno 
however finds a different expression for the deflection. Moreover he computes how the 
time delay is affected by the cosmological constant, which is particularly interesting in 
the light of the impressive progress [5] in the observation of such a time delay. 

The aim of the present paper is to to give a detailed, self-contained account of how 
to compute the deflection in presence of a cosmological constant and only in terms of 
quantities observed in lensing clusters. In particular we do not use the lens equation. Our 
computations are to first order in Newton's constant and in the deflection angle. We do 
not neglect higher terms in the cosmological constant as they do contribute significantly 
in certain lensing clusters. 



2 The set up 

Consider a universe, which is empty except for one static, spherical, non-rotating mass 
M, the lens L. A source S, at rest with respect to the lens, emits photons, which are 
observed at nostra terra T also assumed at rest. We neglect the masses of the source 
and of our local group. We use polar coordinates (r, 9, (p) centered at the lens. Because 
of spherical symmetry the photons' trajectory is in a plane that we take to be 9 = tt/2. 
The angle if is measured with respect to the axis defined by the source, <ps = 0. With a 
cosmological constant, the gravitational field outside the mass M is given by the Kottler 
metric, 

dr 2 = Bdt 2 -B- l dr 2 -r 2 d V \ 9 = n/2, B = l-^--\Ar 2 . (1) 

Our initial conditions for the photon are position (r^, (ps = 0) and velocity characterised 
by the coordinate angle es '■— rs\dip/dr (rs)\- Let r be the peri- lens. We will suppose 
2GM/tq <C 1 and e$ <C 1 and keep only terms linear in these two quantities. We will 
also suppose that Ar 2 /3 < 9/10 to avoid the coordinate singularity at the equator of the 
de Sitter sphere. The final conditions are here on earth (tt, <Pt) where the photons are 
detected together with their (coordinate) angle := r T |d<^/dr [tt)\- 

We are interested in situations where more than one images of the source are seen. 
Because of the spherical symmetry there can be only one or two images depending on the 
distances of earth and source from the lens. There is however one exception, when earth, 
lens and source are aligned. Then = r$, v?t — tt and we get an infinite number of 
images forming Einstein's ring. 
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We will concentrate on the non-aligned case and denote by e' s , r' , e' T the parameters 
of the lower trajectory, figure 1. 




x 



Figure 1: A double image 

Our tasks are: 

• Integrate the photons' geodesies and compute <^r as a function of r s , es, t t and 
(p' T as a function of rs, e' s , r?- The surprising, but well known, fact is that in this 
calculation the cosmological constant drops out, which is not the case for massive 
geodesies. 

• Compute the coordinate distance r T as a function of r s , ex, e' T , such that the two 
trajectories meet, ip' T = (fx- 

• Compute the physically measured angles «t and a' T from the coordinate angles €t 
and e' T . Here the cosmological constant re-enters the scene [Tj. 

• Compute the area distance di of the lens as seen from earth from the coordinate 
distance tt- 

• Compute the area distance ds of the source as seen from earth from the coordinate 
distances Tt, rs and the angles «t and a' T . Here the cosmological constant will also 
play a role. 

Finally we will compare the theory to the observation of the lensing cluster SDSS J1004 
+4112 with a quasar as source. The angles ax and a' T are measured as well as the mass 
of the lens and the redshifts z L of the lens and z s of the source [6j [7]. 

• We will compute the area distances di and ds from their redshifts zl and zs using 
the Hubble diagram for a flat 3-space with cosmological constant and 27 % of matter. 



3 Integrating the geodesies 

We start with the list of the non-vanishing Christoffel symbols for the Kottler metric with 
6 = 7r/2 and denote ' := d/dr, 

T t tr = B'/(2B), r tt = BB'/2, T r rr = -B'/(2B), (2) 
T r w = -rB, r% = l/r. (3) 
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The geodesic equations read: 



i+B'/Btr = 0, 
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r + \BB'i 



\BB'r 2 



rBip 2 = 0, 



ip + 2r rip = 0, 



(4) 
(5) 
(6) 



where we denote the affine parameter by p and ' := d/dp. We immediately get three first 
integrals: 

i = l/B, r 2 (p = J, r 2 /B + J 2 /r 2 - 1/B = -E. (7) 

The last two come from invariance of the metric under rotations and time translations and 
the integration constants J and E have the meaning of angular momentum and energy 
per unit of mass. For the photon, E = 0. Eliminating affine parameter and (coordinate) 
time we get: 

dr 

dip 

At the peri-lens, dr/dip(r ) = and therefore J = r -B(r ) -1 / 2 . Substituting J into 
equation (jSJ), the cosmological constant drops out and we have: 

'1/2 



±r^r 2 /J 2 - B. 



= ± 

dr 



r \7 r2 / r c 



2GM 2GM 



r r + r 



(9) 



From now on we will omit terms of order (GM/ro) 2 , e| and (GM/ro)es and write equalities 
up to this order with a ~ sign. For example: 

ts = r s \dvldr{r s )\~(r 2 s /rl-l)- 1 f 2 



ra/r s . 



(10) 



Note that for the upper trajectory, dip/dr is negative for r between rg and r , positive 
between r and r T . Therefore 



ip T 





dip 


dr + 




dip 


dr. 


fro 


dr 


Jro 


dr 





fill 



Using J x 1 (x 2 — 1) 1 I 2 dx = — arcsinl/a;, J x 2 (x 2 — 1) x l 2 dx = (x 2 — \) x l 2 jx, J(x + 
l) _1 (.x 2 — l) _1 / 2 d.x = [(x — l)/(x+ l)] 1 ^ 2 , we get to linear order: 

/ r s \ AGM , . 

ip T ~ 7T - e s 1 + — + . 12 

V r Tj e s r s 

For the lower trajectory the signs of dip/dr are opposite and we have 

AGM 



ip' T ~ 7r + e' s 



1+2 

r T 



(13) 



We can trade es for ex using r ~ es r s ~ etTt and likewise for the primed quantities. 
For the two trajectories to meet, ip? = ip' T , we must have: 

r T AGM 
_Z i. 

r s e T e' T r T 



(14) 
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4 From coordinate angles to physical angles 



Let us compute the relation between the coordinate angle ex and the corresponding phys- 
ical angle a,?. Officially since October 1983, lengths are measured in terms of proper 
time of flight of photons, one nano-second ~ 0.3 m ~ 1 foot. The coordinate time 
dt it takes the photon to travel from (rx, (px) to (rx — dr, (px) is computed from = 
B{r T ) dt 2 — 5(r T ) _1 dr 2 , see figure 2. 




Figure 2: Measuring the angle «t in ns/ns 
Likewise the coordinate time dt it takes the photon to travel from (rx — dr, (px) to 

~ 2 

(rx — dr, (fx — dip) is computed from = B(rx) dt — r T dip 2 . The we have, 



dr dtWB(rx) rxdip , — - — r / „ 

tan«T = ~r~ = — ^=^= = ; =fTV Birx) ~ £r\/ 1 — Ar4/3 ~ ctx- (15) 

dr dt^B^x) dr/^W^x) V 



Similarly we get a' T ^ e' T yl — Ar|,/3. Rindler & Ishak [lj compute the physical angle by 
means of the usual formula 

= (dr,0)-(dr,dy) 

C ° SC * T [(dr,0) • (d^O)] 1 ^^^^) ■ (dr,d^)]V2 ' 1 > 

where the scalar product ■ comes from the negative of the spatial part of the Kottler 
metric. The two ways to compute the physical angle agree to all orders. 

Finally we note that in Schwarzschild's solution and far out in the asymptotic region, 
physical and coordinate angles coincide. 



5 From coordinate distances to area distances 

Imagine a standard candle radiating photons isotropically at the position of the lens. Its 
area distance ^ as seen from the earth is defined by the relation between the infinitesimal 
solid angle dVt in which photons are radiated and the infinitesimal area dS* of the light- 
sensitive plate on which they are collected on earth, dS/dfl = And 2 . The isotropy of the 
standard candle is of course with respect to physically measured angles a which in this 
position do coincide with the polar angles 6 and <p, see figure 3. 




L 



Figure 3: Area distance of the lens as seen from earth 
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The area is again measured by means of the time of flight of test photons traveling along 
its edges. For the edge in the plane 9 = 7r/2 this coordinate time of flight is given by 
= B(r T )dt 2 -r\ dp 2 . The proper time is therefore dr = r^dip. For symmetry reasons 
we get the same proper time for an edge orthogonal to the plane 6 = n/2 and consequently 

d L = r T . (17) 

For the area distance d$ of the source as seen from the earth, the calculation is more 
involved. The isotropy of the standard candle in the plane 9 = n/2 is defined with respect 
to the physical angle as whose period is 2ti B{rs) 1 ^ 2 • The proper time of flight length of 
the infinitesimal edge at earth in this plane is d£\\ = rr|d<^r|, see figure 4. 




Figure 4: Area distance of the source as seen from earth 
Differentiating equation (fl2|) with respect to es and using equation (Tl~4"I) we get 

ilpT „ -( 1+ r JL)( 1+ 5:) des ^_( 1+ rjL)( 1+ <\ ^ 



r T J \ e T ) \ r T J \ a T J ^jB(r s ) 



r T J V oi T J y/\ - Arf/3 

For an edge orthogonal to the plane 6 = tt/2, the calculation is easier because of the axial 
symmetry around the axis lens - source. An infinitesimal rotation by a coordinate angle 
dr] of period 2ir results in an infinitesimal length at the earth of 

( r T \ \otT — C^t\ 
1+ — — dr]. (19) 

r s J y/B(r T ) 

Replacing the coordinate angle dr? by a physical angle we finally get the area distance 
from dS = d£\\d£ i : 



y/l - Ar 2 /3 y V a T ) r s ^1 - Ar^/3 

Note that a singularity occurs when earth, lens and source are aligned, olt = ot' T . This 
comes from the focusing of the lens. 
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6 From redshifts to area distances 



Here comes the shaky part of the reasoning. Indeed, as long as we do not have a solution of 
Einstein's equation interpolating between Kottler's and Friedmann's solutions, we do not 
know above what length scale the masses of the other galaxies and their expansion must 
be taken into account. Justified only by this ignorance, we make the crude assumption 
that the measured redshifts are exclusively due to expansion, although we had put the 
source and the earth at rest in the Kottler metric. We will use the Hubble diagram 
from the standard model of cosmology with a flat 3-space, a cosmological constant A = 
1.5 • 10 -52 m~ 2 , a Hubble parameter Hq = 2.5 • 10~ 18 s _1 and 27 % of matter in order to 
compute the area distances. A Runge-Kutta integration of Einstein's equation is used to 
obtain the scale factor a(t). For numerical convenience, we choose the initial condition 
a = 1.2 ■ 10 26 m. Then the area distance d follows from the redshift z: 

f z da(t(z))/dt .„ 
d(z) = / \ /J' dz, 21 
Jo a(t{z)) 1 ; 

where by abuse of notation we write t(z) for the inverse function of z(t) := l/a(t) — 1. A 
rough fit to the numerical solution is given by [8], 

a{t) ~ a (pH t) 1/p , p = 0.69 . (22) 

It is good to 3 % up to z = 2 and yields 

[l-p)H Q 

Because of its singularity we cannot use the area distance ds (120]) to estimate the 
position of the source. The singularity is interpreted as magnification rather than an 
actual get-together. But then we need another assumption to estimate the position of the 
source. The simplest assumption coming to mind is 

ds = r + rS (24) 
A/1 " ArI/3 

which would be the correct area distance in absence of the lens. Of course one can argue 
that when we send the mass of the source to zero, we should send A to zero at the same 
time, d s = r T + r s . 



7 SDSS J1004+4112 

Let us see how our assumptions compare to observation. Consider the lensing cluster of 
SDSS J1004+4112 and the quasar as source [HI [7j. As we have at least 4 images, the 
cluster cannot be spherically symmetric. We will again close our eyes and consider only 
the images C and D with ay = 10" ± 10 % and a' T = 5" ± 10 %. The mass of the cluster is 
M = (1±0.2)-10 44 kg. The cluster has a redshift of Zl = 0.68 yielding d L = r T = 7.0-10 25 
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m from the numerical integration. For the quasar we have zs = 1.734 and ds = 13.7- 10 25 
m, which translates into rs = 7.5 • 10 20 m with the singular area distance equation (1201) . 
However with this value the condition (I14p 

r T AGM 4 9 , . 

— ; — 1 - Ar£/3) - 1. 25 

r s a T oi T r T 

for the two photon trajectories to meet cannot be satisfied whatever the value of A, zero 
or positive. If we use the non-singular area distance (l24"j) we get rs = 5.6 ■ 10 25 m and 
the condition for the trajectories to intersect can be met within the experimental error 
bars for the cluster mass M and the angles ay and a' T only if the cosmological constant 
satisfies 

A > 0.81 • HT 52 m" 2 . (26) 

On the other hand, if we assume ds = ty + rs we get rs = 6.7 ■ 10 25 m and the condition 
for the trajectories to intersect can be met within the error bars only if the cosmological 
constant satisfies 

A = (2.5 ± 1.5) • KT 52 m" 2 . (27) 

In both cases we have not varied the cosmological constant in the Hubble diagram and 
it is encouraging that both constraints, (126]) and (1271) are compatible with the present 
observational bounds, A = (1.5 ±0.7) ■ 10 -52 m~ 2 , from the standard model of cosmology. 
Note that Ishak, Rindler et al. |2], presumingly assuming ds = Tt + rs, find an upper 
bound from two other lensing clusters. 



8 Conclusions 

We agree with Rindler & Ishak, gravitational lensing depends on the cosmological con- 
stant. Our formula (1251) for this dependence agrees with theirs |T] in the aligned case. We 
also agree with them j2] , this dependence is not negligible for certain clusters and further 
confrontation with observational data is necessary. 

Of course the main problem remains to find a solution of Einstein's equation that 
interpolates between a homogeneous family of static, curved Kottler solutions and the 
expanding, flat Friedmann solution. We hope that Rindler & Ishak's ground breaking 
work will give new impetus to this fundamental problem. 

Many years ago, dark energy was proposed as an alternative to the cosmological con- 
stant. And still, nobody has ever told us how dark energy modifies the Schwarzschild 
solution. 

Acknowledgements: It is a pleasure to thank Christoph Stephan for lively discussions. 
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